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Abstract 



We investigate the independent points in the sky require to search the periodic gravita- 
(2JT)! tional wave, assuming the noise power spectral density to be flat. We have made an analysis 

with different initial azimuth of the Earth for a week data set. The analysis shows significant 
difference in the independent points in the sky for the search. We numerically obtain an ap- 
proximate relation to make trade-off between computational cost and sensitivities. We also 
discuss the feasibility of the coherent search in small frequency band in reference to advanced 
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1 Introduction 

The first generation of laser interferometric gravitational wave observatory (LIGO) [Abramovici 
(1992)] is in operation. However, its has been planned that by the end of the initial LIGO science 
run it would undergo upgrades to significantly improve the sensitivity. Now referred to as advanced 
LIGO (formerly LIGO II) [Weinstein, 2002]. The detector will sweep their broad quadrupolar 
beam pattern across the sky as the earth moves. Hence the data analysis system will have to carry 
out all sky searches for its sources. In this, search of the PGW without a priori knowledge appears 
to computationally quite demanding even by the standard computers expected to available in the 
near future. It appears that due to limited computational resource it will be not feasible to perform 
all sky all frequency search in the months/year data set. However, if advanced LIGO achieve 
its design sensitivity ~ 10~ 23 or better (Weinstein, 2002), then it may be feasible to perform all 
sky search for a day to week data set in small frequency band for the sources emitting signal 
of amplitude > 1CT 26 . The search of the potential sources may be more significant, if done in 
the frequency band where most of the Pulsars are detected by other means. Also, the choice of 
sophisticated, optimal data analysis methods and a clever programming is also integral part to 
search the signal buried in the noise with the available computation power. 
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The current status of the search indicates that its important to detect the gravitational waves 
(GW) rather finding the source location more accurately. Hence, one would like to do minimum 
Doppler correction or/and to make the search templates. In reference to the all sky search, Schutz 
(1991), has introduced the concept of patch in the sky as the region of space throughout which 
the required Doppler correction remains the same. He roughly estimated the number of patches 
required and shown that for 10 7 sec. observation data set of one KHz signal would be about 1.3 x 
10 13 . This also indicates the bank of templates require for the coherent all sky search. The size of 
the patch may be increased, hence reducing the number of points require for the Doppler correction, 
by manipulating the output of the detector. Which in turn demand the detail investigation of the 
parameters affecting the phase of the modulated signal. In this, the initial azimuth of the Earth play 
a vital role in the modulation of the signal, particularly for the analysis of day to week data set. 
Hence, in the next section, we incorporate the initial azimuth of the Earth in the Fourier transform 
(FT) obtained by Srivastava and Sahay (2002a,b) for arbitrary observation time. 

The limited computation power make to search the signal in the short observation data set, may 
be a day /week, so in section 3, employing the concept of fitting factor (Apostolatos, 1995), we 
investigate the independent points in the sky (N sky ) require for an all sky search for a week data 
set with different initial azimuth of the Earth, assuming the noise power spectral density to be flat. 
Also, we numerically obtain an approximate relation to make trade-off between computational 
cost and sensitivities. In section 4, we discuss the feasibility of the coherent search in reference to 
advanced LIGO. Section 5 contains the conclusions of the paper. 

2 Fourier transform 

The FT analysis of the frequency modulated PGW signal has been done by Srivastava and Sahay 
(2002a,b) by taking account the effects arising due to the rotational as well as orbital motion of 
the Earth. However, they have neglected an important parameter, the initial azimuth of the Earth, 
which affects significantly in the spacing of the parameters space for an all sky search. To obtain 
the FT with taking account the Earth initial azimuth, we rewrite the phase of the received PGW 
signal of frequency f a at time t given by them and may be written as 

$0) = 2nf t + Z cos(a£ rot — a) + A/"cos(£ rot — 5) — M (1) 
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where 0, 0, R e , R se , w orb w rot , a and e represent respectively the celestial co-latitude, longitude, 
Earth radius, average distance between Earth centre from the origin of SSB frame, orbital and ro- 
tational angular velocity of the Earth, co-latitude of the detector and obliquity of the ecliptic. Here 
(3 or b and f3 rot are the initial azimuth of the Earth and detector respectively. 

To estimate N sky , it is sufficient to consider either of the two polarisation of the signal given as, 



hence, we consider the '+' polarisation of amplitude unity, given as 



where h 0+ , h 0x are constant amplitude of the two polarizations. 

Now let us assume h(t) to be given on finite time interval [0, T obs ] assumed to be the observation 
period. Now its straight forward to obtain FT in the similar way as obtain by Srivastava and Sahay 
(2002b), and may be given as 



h+(t) = h 0+ cos[$(t)] 



(3) 



h x (t) = h 0x sin[$(t)] 



(4) 



h(t) = cos[$(f)] 



(5) 




(6) 
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where J stands for the Bessel function of first kind. Using the symmetrical property of Bessel 
function we reduce the computation time appreciably by rewriting h(f) as 
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(8) 



X = sin(M - mix j2) 

y = cos(M - mn/2) 

U = sin z/£ D cos m(£ D — 5) — ™ {cos z/£ D sin m(£ — 5) — sin m5} 

V = cos z/£ cos m(£ Q — 5) + 77 sin z/£ Q sin m(£ — 5) — cos m<5 

The transform contains double infinite series of Bessel function. However, for analysis the 
order of Bessel function required to compute h(f) in the infinite series are given as (Srivastava and 
Sahay, 2002b). 

k « 3133.22 x 10 3 sin 9 ( J , (10) 

\lkHz J 

The accuracy and range of validity for large order and/or argument has been discussed by 
Chishtie et.al (2005). 
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3 Independent points in the sky 



The study of the independent points for an all sky search has been made by many research workers 
[Schutz (1991), Brady et. al. (1998), Brady and Creighton (2000), Jaranowski and Krolak (1999, 
2001), Astone et.al (2002)] for the coherent and/or incoherent search. The coherent search means 
cross correlating the data with the bank of search templates. While incoherent search implies 
adding of the power spectra by dividing the data into N subsets, performing a full search for each 
subset, and adds up the power spectra of the resulting searches. In this case, there is loss in S/N 
ratio by a factor of \fN in relation to coherent search as power spectra are added incoherently. 
However, irrespective of the method of search, the optimal spacing in the (9, 4>) parameters for an 
all sky search is a problem of interest. 

For the coherent search, one have to make bank of search templates to detect the signal. The 
bank of search templates are discrete set of signals from among the continuum of possible signals. 
Consequently all the signals will not get detected with equal efficiency. However, it is possible to 
choose judiciously the set of templates so that all the signals of a given amplitude are detected with 
a given minimum detection loss. Fitting factor (FF) is one of the standard measure for deciding 
what class of wave form is good enough and quantitatively describes the closeness of the true 
signals to the template manifold in terms of the reduction of S/N arising due to the cross correlation 
of a signal outside the manifold with the best matching templates lying inside the manifold, given 
as 



FF = W)\hr(f;6 T ,<h)) u 

V(M/; Or, <h)\hr(f\ Br, <h))W)\h(f)) 

where h(f) and hr(f; 9t, 4>t) represent the FTs of the actual signal and the templates respectively. 
The inner product of two signal hi and h 2 is defined as 

- 2 / _ df 
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where * denotes complex conjugation and S n (f) is the spectral noise density of the detector. 

To estimate optimal spacing in the parameters space one have to careful investigate the param- 
eters contain in the phase of the modulated signal. Hence we check the effect for different (3 or b 
in h(f) for a week data set for the LIGO detector at Livingston [the position and orientation of 
the detector can be found in Allen (1995)] of unit amplitude signal for f = 50 Hz and (9, 4>) 
=(7r/18, 7r/4). We take the ranges of k and m as 1 to 27300 and -15 to 15 respectively. We found 
that Earth azimuth affects the FM spectrum severely. Hence, for the coherent search, we investigate 
its effect for the number of independent of points in the sky. 

To estimate N sky , we consider the LIGO detector at Livingston, receive a PGW signal of fre- 
quency f = 50 Hz from a source located at (9, 4>) = (1°, 45°). First, we chosen the data set such 
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that [3 orb = at t = 0. In this case we take the ranges of k and m as 1 to 2800 and —15 to 15 
respectively and bandwidth equal to 50 ± 3.28 x 10~ 4 Hz for the integration. Now, we select the 
spacing AO = 0.45 x 10 -4 , thereafter we maximize over <p by introducing spacing Acp in the so 
obtained N s k y and determine the resulting FF. In similar manner we obtain N s k y for p or b = n/6 
and 7r/2. The results obtained are shown in the Fig. (Q~|). We also plot the templates spacing A</> 
in the </>— parameter with FF shown in the Fig. ©. Interestingly, the nature of these curves are 
similar, we have obtained a best fit of the graphs, given as 

N sky = 10 18 [a o + a x x - a 2 x 2 + a 3 x 3 - a 4 x 4 + a 5 x 5 - a 6 x 6 + a 7 x 7 }; 0.85 < x < 0.995; (14) 

FF = b Q + biy - b 2 y 2 + hy 3 ; 0.037 < y < 0.69; (15) 
where and b Q ..b 3 are constants as given in Table 1 and 2 respectively. 
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Table 1: 


Coefficients of the best fit graphs obtained for the N s 


ky with FF. 






Porb 




61 




h 













0.994123 


0.149451 


3.94494 


2.64819 










7r/6 
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2914.21 


37817.3 







Table 2: Coefficients of the best fit graphs obtained for the FF with A</>. 



In view of the above investigation, the spacing A0 in the 0— parameter may be expressed as 

A0 = G(FF, f , 6, 0, Tobs, Porb) (16) 

Equation (TT5T ) is a third order polynomial, hence complicates the solution. However, one would 
like to do data analysis for FF > 0.90. Therefore from the figure © we obtain a very good 
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dependence of FF for minimum N sky and may be given as 

FF = -2.85657 A <p 2 + 0.00230059 A + 1.00018 (17) 
From Eqs. ([Tot , and (TT71 ), we may write 

T{FF, 50, 1°, 45°, lu;, /5 orfe ) « 4.02684 x 10" 4 ± 0.591667^1-00018 - FF (18) 

Above equation can be relevant to make trade-off between computational costs and sensitivities 
i.e. for the selected FF one can estimate N sky . However, there is no unique choice for it. Here 
we are interested in the estimation of N sky such that the spacing is maximum resulting into the 
least number of points require for the search. As mentioned earlier, there is stringent requirement 
on reducing computational costs. Accordingly, there is serious need of adopting some proce- 
dure/formalism to achieve this. For example, one may adopt the method of hierarchical search. 
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Figure 1 : N sky with FF at different (3^. 
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Figure 2: FF with A0 at different [3 orb . 

4 Computational costs 

In view of the above analysis it will be interesting to know the feasibility of the all sky search with 
the target sensitivity of the advanced LIGO. The computational costs of the data analysis basi- 
cally depends on the floating point operations (flops) require to perform the Fast Fourier transform 
(FFT). Hence in terms of FFT, the flops for the data reduction upto frequency / for T obs of the 
interferometer output may be given as (Press et al., 1986) 



N flops 

2fT obs log 2 {2 fmaxTobs ) 



(19) 



Now for the given mismatch (FF), if N p is the number of independent points to perform all sky 
search, then the flops will be 



N flops (FF,N p ) = 2fT obs N p log 2 (2f max T obs ) 



(20) 



Hence, for the mismatch of 3% and without manipulating the data in reference to f3 orb , the flops 
for the search of PGW signal upto 50 HZ in the week data set will be 2.38 x 10 19 , 2.58 x 10 20 and 
6.23 x 10 20 for f3 orb = 0, 7r/6 and n/2 respectively, assuming other operation need negligible flops 
compare to FFT. However, the lower cut off frequency of the LIGO I/II is 10/40 Hz. Hence the 
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analysis shall be done above the lower cut off. Also, the search will be more significant if one 
perform in the most sensitive band of the detector. In this, if one would like to perform all sky 
search in a small band say 5 Hz then the minimum flops for the on-line analysis (a week data gets 
analyzed in a ~ week time) will be 3.44 x 10 12 . The flops require may be further reduce, if one 
perform hierarchal search. Hence, it may be feasible to perform limited frequency all sky search 
of signal amplitude > 1CT 26 in the output of such a sensitive detector with a ~Tflops computer. 



5 Conclusions 

In this paper we have incorporated the initial azimuth of the Earth in the FT of the frequency 
modulated PGW signal and investigated its effect in the independent points in the sky require for 
the search of PGW. We found that the N sky for the search in the output of one week data set 
varies significantly with (3 orb . For the case investigated here, we observe that for FF = 0.97 
approximately 1.53 x 10 10 , 1.66 x 10 11 and 4.0 x 10 11 N sky will be require when f3 orb = 0, n/Q 
and 7r/2 respectively. Hence, the analysis may be useful to reduce the computational cost for a 
coherent all sky search. However, also the inspection of the phase of the modulated signal reveals 
that reduction in N sky depends on time scale of integration, shorter the T obs more the difference in 

The reduction in N sky is large, so we studied the feasibility of all sky search in reference to 
advanced LIGO and found that in the band of 5 Hz one may perform on-line all sky search of the 
PGW signal of amplitude > 10~ 26 with a ~ Tflops computer. The relation given by equation (fl~8l) 
may be useful to make trade-off between computational costs and sensitivities for the search of 
periodic gravitational waves. The issue to reduce the flops for the all sky search is a problem of 
interests and hence need more studies/investigations. 
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